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Abstract 

We construct Q-operators for the open spin-^ XXX Heisenberg spin chain with 
diagonal boundary matrices. The Q-operators are dehned as traces over an 
inhnite-dimensional anxiliary space involving novel types of reflection oper¬ 
ators derived from the bonndary Yang-Baxter eqnation. We argne that the 
Q-operators dehned in this way are polynomials in the spectral parameter 
and show that they commute with transfer matrix. Finally, we prove that the 
Q-operators satisfy Baxter’s TQ-equation and derive the explicit form of their 
eigenvalues in terms of the Bethe roots. 
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1 Introduction 

Integrable spin chains are prominent examples of quantum integrable models. They bear a 
close relation to two-dimensional integrable quantum field theories as well as lattice models 
and rest on surprisingly rich mathematical structures. Furthermore, integrable spin chains 
play an important role in condensed matter physics and in the AdS/CFT correspondence. 

The key example of an integrable spin chain is the periodic spin-| XXX Heisenberg spin 
chain. This one-dimensional quantum system describes the interaction of neighbouring 
spins aligned on a line with closed periodic boundary conditions and has been solved by 
Hans Bethe pQ. 
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Curiously, the spin chain need not have closed boundary conditions to be integrable. 
In fact, it was shown in |2j that the nearest-neighbour Hamiltonian for the Heisenberg 
spin chain with open boundary conditions 

L-l 

H = ^ (Ji ■ (Ti+i ^ , (1.1) 

i=l 


can be solved using the coordinate Bethe ansatz. Here ^ and ^ are arbitrary boundary 
parameters, L denotes the length of the spin chain and at the Pauli matrices at spin chain 
site i acting on the 2^-dimensional quantum space. 

The quantum inverse scattering method is a powerful tool to study integrable models. 
It allows the transfer matrix to be constructed using R-matrices that are solutions of 
the Yang-Baxter equation. An introduction, as well as references, can be found in [3]. 
Initially formulated for closed systems, the quantum inverse scattering method was further 
extended to integrable models with open boundaries by Sklyanin [1]. In particular, this 
allowed the transfer matrix to be constructed for the open XXX Heisenberg spin chain. As 
for the closed case the method relies on the Yang-Baxter equation but in addition to the 
R-matrices also incorporates certain boundary matrices obeying the boundary Yang-Baxter 
equation. 

Baxter Q-operators are of distinguished importance in the theory of integrable systems. 
They contain the information about the eigenfunctions as well as the Bethe roots. They 
are named after Baxter who introduced a Q-operator to solve the eight-vertex model |5]. 
Nowadays, interest in the Q-operators stems from various branches of theoretical physics. 
Q-operators are of central importance in the ODE/IM correspondence [B] which relates 
quantum integrable models to ordinary differential equations and play a major role in the 
formulation of the fermionic basis [7j. A systematic, representation theoretical approach 
to Q-operators was recently studied in [HI E]. Furthermore, the results obtained for the 
spectral problem of planar A/" = 4 super Yang-Mills theory, see [TU] and references therein, 
suggest that Q-operators play a central role in the study of the integrable structures 
appearing in the context of the AdS/CFT correspondence, see e.g. HB and references 
therein and thereof. 


Even though the Q-operator is often referred to as the most powerful tool for the exact 
diagonalisation of integrable systems it has only been studied for a few examples. In fact, 
the construction for the periodic XXX Heisenberg chain in the presence of a magnetic field, 
along the lines of ca. is rather recent HBJ. See also m, where these Q-operators have 
been diagonalised using the algebraic Bethe ansatz. Apart from the construction in terms 
of integral operators for certain non-compact open magnets that in particular appeared in 
the context of QCD [13 US] little is known about Q-operators for open spin chains. 

In this article we construct the Q-operators for the open XXX spin-^ Heisenberg spin 
chain. Following the quantum inverse scattering method we construct the Q-operators as 
traces of single-row monodromies similar to the ones that appeared in [T3] and novel types 


^Throughout this article we restrict ourselves to the case of diagonal boundary matrices with the 
Hamiltonian given in (1.11. 
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of boundary operators which are derived from the boundary Yang-Baxter equation. We 
show that the Q-operators commute with the transfer matrix, prove Baxter’s TQ-equation 
on the operatorial level and derive their eigenfunctions. 

The article is organised as follows. We hrst review the construction of the transfer 
matrix for the open XXX spin-^ Heisenberg spin chain in Section In Section we 
define the Q-operators. We introduce the Lax operators for the single-row monodromies 
and derive the boundary operators for the Q-operators from the boundary Yang-Baxter 
equation. Section |4] discusses various properties of the two Q-operators. We show that the 
transfer matrix commutes with the Q-operators and elaborate on the symmetry which 
relates them. Furthermore, we carefully study the trace in the auxiliary space of the 
Q-operators and show that it yields polynomials in the spectral parameter. More precisely, 
we demonstrate that the entries of the Q-operator are arranged in blocks. We evaluate them 
at leading order of the spectral parameter. In Section]^ we prove Baxter’s TQ-equation 
on the operatorial level and argue that the eigenfunctions of the Q-operators yield the 
desired Q-functions. Finally, we discuss the results and provide an outlook in Section 


2 The open Heisenberg spin chain 

In this section we review the construction of the fundamental transfer matrix of the open 
XXX Heisenberg spin chain in the framework of the quantum inverse scattering method 
following |1]. 

At first we define the fundamental R-matrix 


R{z) = z + F, (2.1) 

with the permutation matrix P = Yl\b=i ^ab ® ^ba- Here z is the spectral parameter and 
Cab denotes the 2x2 matrices with the entry “1” in the row and column and “0” 
elsewhere. Thus, the R-matrix acts on the tensor product of two spaces 


R{z) : ® ^ ® 


( 2 , 2 ) 


The R-matrix in (2.1) is a solution to the Yang-Baxter equation 


Ruix - ii)Ri3{^ - 2)i?23(y - 2 ) ^ Ri^iy - y)Ri3(x - z)Ri 2 (x - y ), 


(2,3) 


where we defined 

Ri 2 iz) = R{z) ® I, (2.4) 

and Ri 3 , R 23 , correspondingly. Furthermore, one hnds that the R-matrix satishes the 
unitarity and crossing unitarity relation 


R{z)R{-z) = 1 - R*' {z)R^^ {-z - 2) = -z{z + 2). 


( 2 . 6 ) 


Here ti denotes the transposition in the first space of the symmetric R-matrix, cf. (2.2). 
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In order to construct the transfer matrix of the open spin chain we dehne the quantum 
Lax operators which coincide with the R-matrices introduced above 

C^^{z) = Ro,{z). (2.6) 

They act on the spin chain site in the quantum space and on the two-dimensional 
auxiliary space denoted by “0” which has been suppressed in the notation above. The 
transfer matrix T of the open spin chain is dehned from the double-row monodromy matrix 
U as 

T{z) = iiolC{z)U{z) with U{z) = M.{z)K{z)jCi[z). (2.7) 

The single-row monodromy matrices M.{z) and jCi(z) are constructed from the tensor 
product of the Lax operators with z = 1,..., L in the quantum space and multiplication 
in the auxiliary space as 

M(z) = (z)£[21 ( z ) • ■ ■ (z), M(z) = (z) ■ ■ ■ (^) . ( 2 . 8 ) 

Here L is the length of the spin chain, i.e. the number of sites that build up the 2^- 
dimensional quantum space. The boundary operators /C and /C are diagonal matrices that 
only act non-trivially in the auxiliary space and read 


/C(z) = 


p + z + 1 0 

0 p — z 





(2,9) 


Here, as discussed at the end of this section, the boundary parameters in (1.1) are related to 


p and q via ^ = 1/p and ^ = 1/q. The boundary operators in (2.9) can be determined from 


the boundary Yang-Baxter equation. In particular, the following boundary Yang-Baxter 
equations hold 


R{x — y){U{x) ® I)R{x + y)(I ®U{y)) = (I ®lA{ii))R{x + y){U{x) ® I)R{x — y) , (2.10) 


and 


{l®K,{y))R{—x — y — 2){K,{x) ®l)R{y — x) = R{y — x){K,{x) ®l)R{—x — y — 2){l®lC{y )), 

( 2 . 11 ) 

where in addition to the explicit form of the boundary matrices (2.9) we used the Yang- 


Baxter equation (2.3) to show (2.10) which involves the double-row monodromy lA. 


The transfer matrix satisfies the crossing relation 


= T( 


1 ), 


( 2 . 12 ) 


This relation is shown in Appendix as it is needed in Section For alternative proofs 
see e.g. [niiiH]. As a consequence of the boundary Yang-Baxter equations above and the 


unitarity relations in (2.5) one hnds that the transfer matrix commutes for different values 


of the spectral parameter 


lT(i),Tfe)]=0. 


(2.13) 
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The spectrum of the transfer matrix can be determined using Bethe ansatz methods. The 
eigenvalues of the transfer matrix T are given in terms of Baxter Q-functions by the 
TQ-eqnation 


T{z) 


{z±p){z±q) 


{z + lf^+^ Q±{z-l) 
z+^ Q±(z) 


+ (zTP+l)(zTg + l) 


z^^+^ Q±(z + 1) 
z + ^ Q±(z) 

(2.14) 


The TQ-equation can be obtained from the algebraic Bethe ansatz. Depending on the choice 
of the reference state, all spin up or all spin down, the eigenvalues of the transfer matrix 
are either parametrised by the Q-fnnction or Q-. The Q-fnnctions are polynomials 
of degree 2m± parametrised by the Bethe roots zf. Here m± is the magnon nnmber and 
denotes the nnmber of np (down) spins in the sea of down (np) spins. The Q-fnnctions 
can be written as 


m± 


Q±(z) = 


n(2-y)(2+y+i) 


(2.15) 


2m± - LTPTq f.i 

We note that commonly the Q-functions are dehned without the ^-independent normalisa¬ 
tion which drops from the TQ-equation. However, the Q-operator construction presented 
in the next section naturally incorporates the prefactor above. The Bethe equations follow 
from Baxter’s TQ-equation. They take the form 


(zt ±p)(zt ±q)(zt + 1) 


2L 


(Zi =FP + I)(zi =Fq + l)(zi ) 




m± 

n 

k=l 






+ i)(j: 


± 




+ 2 ) 








(2.16) 


Finally, the nearest-neighbonr Hamiltonian (1.1) is related to the transfer matrix via 
T'(0) = 2pq (L + H) and consequently belongs to the family of commuting operators. Here 
we identihed ^ = 1/p and ^ = 1/q. Thus, the energy eigenvalues can be obtained in terms 
of the Q-functions from the derivative of T at z = 0. 


3 Q-operator construction 

In this section we dehne the Q-operators Q± of the open Heisenberg spin chain. We hrst 
introdnce the Lax operators for the single-row monodromies and discnss some of their 
properties. Then we derive the bonndary operators relevant to constrnct Q-operators 
from the bonndary Yang-Baxter equation and hnally dehne the Q-operators of the open 
Heisenberg spin chain. 


3.1 Degenerate solutions to the Yang-Baxter equation 


As for the construction of the Q-operators for closed spin-4 Heisenberg spin chain 531 we 
introduce the Lax operators 
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(3.1) 










Here the pair of oscillators (a, a) obeys the usual commutation relations and N is the 
number operator 

[a, a] = 1, N = aa. (3.2) 

The Lax operators above are solutions of the Yang-Baxter equation 


R{x - y){L±{x) ® I)(E 0 L±(?/)) = (I (» L±(j/))(L±(a;) ® I)R{x - y ), (3.3) 


where R is the fundamental R-matrix dehned in (2.1). In addition to the Lax operators 
L± in (3.1) we define 


U(z) = 


z-N 

a 


U = 


— 1 a 
a z — N 


(3.4) 


Together with the Lax operators in (3.1) they satisfy the unitarity and crossing-unitarity 
relations 


L±(;^)Ld 


-z = 


2)L 


z = 


1 , 


(3.5) 


c.f. (2.5). Here the superscript t denotes the transposition of the 2x2 matrix which leaves 


the oscillators untouched. Note that the Lax operators L± can be related to Lj. by a 
particle-hole transformation and a subsequent multiplication of a diagonal matrix. As the 
particle hole transformation does not change the commutation relations and the tensor 
product of any two diagonal matrices commutes with the R-matrix, we hnd that the Lax 


operators L± satisfy the same Yang-Baxter equation as L±, cf. (3.3). 


3.2 Boundary Yang-Baxter equation 

In analogy to the transfer matrix construction we define the boundary operators 

r(^p-^-l-N) p±{z) 


KM) = 


KM) = 


(3.6) 


pM) ’ r{±q-z-N)^ 

which act diagonally in the oscillator space. Here, as discussed below, the normalisation 
reads 

P±{z) = T{tp-z), p±(z) = r(±g - x) . (3.7) 

We did not find these boundary operators in the literature. The operators K± and Ki are 
the counterparts of the boundary matrices in (2.9) and can be determined through the 
boundary Yang-Baxter equations 


L±(i - !/)K±(2 :)L±(i + y)K(y) = K(y)L±(x + !/)K±(i)L±(i - y), 


(3.8) 


and 


K±{y - a;)K±(a;)L±(-a; -y- 2)lC{y) = lC{y)K±{-x -y- 2)K±(x)L±(?/ - x). (3.9) 

Using the explicit expressions for boundary matrices /C and /C as well as the Lax operators 
L± and L-t we find that this equation is trivially satisfied for the diagonal entries in the 
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2x2 space. Here we additionally assumed that the boundary operators K± and K± are 
diagonal, i.e. only depend on the number operator N = aa. The remaining two dehning 


relations arising from (3.8) are conjugate to each other. The same holds true for the 


dehning relations obtained from (3.9). They read 


aK±(a;) = (±g — a; — 1 — N)K±(x)a and K±(a;)a = aK±(a;)(=Fp — a; — 1 — N). (3.10) 


These equations are naturally solved by the Gamma function as given in (3.6). The 


functions p± and p± cannot be determined from the boundary Yang-Baxter equations (3.8) 


and (3.9). The normalisation in (3.7), which renders the Q-operators to be polynomials in 


the spectral parameter, is obtained in Section 4.4 


3.3 Definition of the Q-operators 

We will now dehne the Q-operators. First, we introduce the double-row monodromies U± 

U±(z) = M±(z)K±(z)M±(z). (3.11) 


Here the boundary operators K± were dehned in ( 3.6[ ). The monodromies M± and M± are 
built in analogy to (2.8) from the tensor product of the Lax operators in (3.1) acting 


on the site of the spin chain and the ordinary product in the oscillator space as 

M±(z) = lW(^)l|'(^) ■ ■ • l.fiz ), M±(^) = ■ ■ ■ lL''(z) . (3.12) 

The double-row monodromy matrix ( |3.11 ) and the boundary operator K± in ( |3.6[ ) allow 
us to dehne the Q-operators Q± as 

Q±( 2 :) = tr K±( 2 :)U±( 2 ;) . (3.13) 

Here the trace is taken in the oscillator space. For convenience we give the explicit form of 


the Q-operators in (3.13). They read 
F(—N — p — z — 1) 


Q+(^) = tr 


T{-p-z) 
F(g - ^) 


a 

^ -t-1 + N 


z -Fl -h N 


and 


Q-(z) = tr 



a 

1 + N 

a 

z-h 1 + N 


(3.14) 


r(-N-?-z) 



(3.16) 


As shown in the next section, these operators are well-dehned and in particular the 
trace yields a hnite result. The entries of the Q-operators are polynomials in the spectral 
parameter arranged in blocks labelled by the magnon number. As also shown in the 
following, the Q-operators commute with the transfer matrix and satisfy Baxter’s TQ- 
equation. 
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4 Properties of the Q-operators 


In the following we discuss some properties of the Q-operators dehned in Section First 
we show that the Q-operators commute with the transfer matrix and discuss a symmetry 
among the Q-operators which is apparent in our construction. Furthermore, we show 
that the trace over the oscillator space is well-dehned and argue that the entries of the 
Q-operators are arranged in L-|- 1 blocks corresponding to the magnon sectors m±. Finally, 
we study the behaviour of the Q-operators when the spectral parameter z is large. We 
hnd that each block is of maximal degree in the spectral parameter and determine 


the overall prefactor apparent in our construction, cf. (2.15). 


4.1 Commutativity 



Figure 1: Schematic proof of the commutativity of T and Q±. 


In this section we show that the Q-operators in (3.13) commute with the transfer 
matrix (2.7). In particular, this is a consequence of the boundary Yang-Baxter equation 
for the double-row monodromies (2.10) and 

L±(a: - 7 /)U±(x)L±(x + y)U{y) = U{y)L±{x + 2 /)U±(a:)L±(x - y ), (4.1) 


as dehned in (2.7) and (3.11) respectively. The boundary Yang-Baxter equation (4.1) 
yields the fundamental commutation relations among the double-row monodromies of the 
Q-operators and the transfer matrix. It follows from the boundary Yang-Baxter equation 
in (3.8) and the Yang-Baxter relations for the single-row monodromies 

AI( 2 /)M±(x)L±(x -y) = L±{x- y)M±{x)M{y ), (4.2) 

and 

A4(2/)M±(x)L±(x - I/) = L±(x - y)M±{x)M{y). (4.3) 

As a consequence of the unitarity relation for the R-matrix ( 2.5[ ) which yields j(4{z) = 
(1 — the RTT-relations (4.2) and (4.3) imply that 


and 


M±(x)L±(x -h y)M{y) = M{y)'L±{x + 2 /)M±(x), 
M±(x)L±(x -h y)M{y) = M{y)L±{x + y)M±{x). 


(4.4) 

(4.5) 
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Now, it is straightforward to show that the Q-operators commute with the transfer 
matrix. Using the explicit form of the transfer matrix and the Q-operators we find 

Cl±{x)T{y) = trotr \]±{x)l{{y)}C{y)K±{x 

trotr [u±(a;)L±(a; y)l{{y)IC{y)L±{-x -y- 2)K±{x) 


X + y + 1 


{x + y + l){x - y) 


trotr [L±(a; - y)\J±{x)L±{x + y)U{y) 

X lC{y)L±{-x -y- 2)K±{x)t±{y - x) 
trotr [l({y)L±{x + |/)U±(a;)L±(a; - y) 

X L±{y - a;)K±(a;)L±(-x -y- 2)}Ca{y) 


(4.6) 


{x + y + l){x - y) 

trotr \u{y)'L±{x + 2 /)U±(x)K±(a:)L±(-x -y- 2)K:{y) 


X + y + 1 
= trotr \u{y)\]±{x)K±{x)K:{y) 

= T(2/)Q±(x). 

Here, we inserted the identity operator in the first and second step via the unitarity and 


crossing unitarity relation (3.5) using the cyclicity of the traces. Employing the boundary 


Yang-Baxter equation for the double-row monodromies (4.1), the boundary Yang-Baxter 


equation in (3.9) and (2.11) for K± and /C and subsequently using again the unitarity and 


crossing unitarity relation (3.5) yields the desired commutation relations 


|T(x),Q±(!/)] = 0 


(4.7) 


The manipulations explained above are schematically depicted in Figure Here each two 
crossing lines denote an R-matrix. The gray vertical lines denote the quantum space while 
the dashed and the black line denote the auxiliary space of T and Q respectively with the 
corresponding K-matrices at the boundaries which are indicated by the black vertical lines 
to the left and right. The hrst and last equality make use of the unitarity and crossing 
unitarity relations as well as the RTT-relations above while the middle equation uses the 
boundary Yang-Baxter equations. 

4.2 Spin-flip symmetry 

For the Heisenberg spin chain the two Q-operators Q± can be related to each other by 
a similarity transformation in the quantum space and a change of sign in the boundary 
parameters p and q. This follows from a relation among the Lax operators L± in (3.1) and 


the explicit form of the boundary operators in (3.6). The Lax operators are related via 



— Lzp(2;). 
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(4.8) 






























Obviously, a similar relation holds for the single-row monodromies M± and M± of the 
Q-operators involving the similarity transformation 


5 = 5 -^ 


0 1 
1 0 



(4.9) 


As a consequence of the symmetry among the single-row monodromies and the form of 
the boundary operators K± and K± in (3.6) one immediately obtains a relation among 
the Q-operators Q±. It reads 




Q=f(^) 


p^-p ■ 

q^-q 


(4.10) 


Finally, we note that the transfer matrix T is invariant under the similarity transformation 
and a subsequent change of sign in the boundary parameters p and q. 


In what follows we often focus on the Q-operator Q+ as it can be related to the Q_ 
using the relation in (4.10). 


4.3 Block structure 


The Hamiltonian as well as the transfer matrix of the open Heisenberg spin chain are block 
diagonal, i.e. they do not contain entries that mix between different magnon numbers. 


This can be seen from the expansion of the boundary Yang-Baxter equation (2.10) at 


order after building the transfer matrix from hl{y), see [3] where the calculation is 

done for the closed chain. We expect that also the Q-operators do not mix between states 
of different magnon numbers such that [Q±(2:),I]^i = 0. This is indeed true. It can 

be verihed by evaluating the boundary Yang-Baxter equation (4.1) at order 


subsequently tracing over the oscillator space after multiplying by K±. However, instead 


of using the fundamental commutation relations that arise from (4.1) we demonstrate this 


property by evaluating the corresponding matrix elements of the Q-operators explicitly. 


This has the advantage of introducing some notation that will be needed in Section 4.5 


From the definition of the Q-operators in (3.13) we find that the matrix elements read 


/ aia2-..CLL 


= tr[K±(^)(M±(^))'^^'"-""k±(^)(M±(^)) 

V / aia2...aL V / 

{c}=l 


b\b2...bL 

Cl C2...C]^ 


(4.11) 


where 


Q±W= E (Q±( 


bib2...bL 

) ^aibi ® ^^2^2 


aia2...aL 


..(8)e„^b^. (4.12) 

Here we dehned the components of the single-row monodromies M± and M± in analogy 


to (4.12). They are given in terms of the components of the Lax operators L-i-(^) = 


J2a,b=i (L±(2^)] ^Gab In the quantum space as 

, bLbi^_i...bi 


(m±(z)) 


ab 

^\bib2...bL 
a\a2...aL 


(m±(^ 






0262 
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In order to show that the Q-operators are arranged in blocks as well as for the discussion 
in Section 14.51 it is convenient to introduce four sets 


= = ^ A k=j} , (4.14) 

of cardinality and i,j G {1,2}. These four sets are in one-to-one correspondence 


to the indices of the quantum space, cf. (4.12), and thus uniquely label a given matrix 


element. From (4.13) and the explicit form of the Lax operators (3.1) we know that the 

monodromy (m+(z)) ^ ^ is a sequence of 1 / 12 ’'^^I operators a, operators a and 

V / 0102...az. 

, bib2...bL 


T{a,c)< 

'22 I 


operators (z -f 1 -|- N). Similarly, (m±(z)) ^ ^ ^ is a sequence of operators a, 

V /ciC2...C£ 


l-^ 2 i’^^l operators a and 1 / 22 ^^I operators ( 2 ; -|- 1 -|- N). As the boundary operators K_|_ and 


K_|_ only depend on the number operator N we hnd that the trace in (4.11) is non-vanishing 
if 

(4.15) 


I r(“-c)| 

|.'l2 I 


I I r{c,'))| _ I r(a,c)| I I r{c,6)| 

“T |.'l2 I — |.'21 I W |.'21 I 


i.e. the total number of creation operators coincides with the total number of annihilation 
M±( 2 ;) ) and (M±( 2 ;)) . The cardinality of the sets can be related 

/C1C2...CZ, V /ClC 2 ...Ci 

to the magnon numbers of the in and out state, i.e. the number of indices with value “2” 
in (ui, 02 ,..., a^,) and (fei, & 2 , • • •, ^l) respectively, via 


_ |7'(“ic)| I |T{a,c)| 

— P 21 I “T r22 I 


b _ I r(c,6)| 


ml = \l 


12 


we find that 


= m+ 


7-(c.b)| 

J-22 \ ■ 

(4.16) 

I 1 r{a.c)| 
“T ^22 1 

— 1 1 1 rFF)! 

~ |.'21 1 “T ^22 1 


(4.17) 


A similar argument holds true for Q_. 

Thus, we find that the Q-operators Q± are organised in L -|- 1 blocks with m± = 
0 ,1 ,..., L where the size of each block is x • 


4.4 Tracing over the auxiliary space 

In this section we show that the trace over the non-vanishing entries of the Q-operators 


is well-dehned and determine the normalisation p± and /5±, cf. (3.7), which renders all 
entries polynomials in the spectral parameter. 

First, we note that the entries of the Q-operators can be expanded in terms of the 
oscillators. Thus, after normal ordering of a given sequence in the monodormy M_|_ and 
anti-normal ordering the sequence in the monodormy M+ as well as using the linearity of 
the trace we hnd that by construction all entries that contribute to the Q-operators in 


(3.13) can be brought to the form 


= tr 

i’i,r2,si,S2 


K±{z) K±{z) a 




(4.18) 
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with ri, r 2 , Si, S 2 G Here, as discussed above, the trace is non-vanishing if there are as 
many creation as annihilation operators 


ri + Si = r2 + S2, 


(4.19) 


cf. (4.15). After inserting the explicit form of the boundary operators (3.6) and reordering 
the oscillators using the commutation relations in (3.2) the argument in the trace of (4.18) 
can be written in terms of the number operator N as 


ri+si 


V, 


ri,r2,si,S2 


H *S'i(ri-1-Si, fc) tr 


r(si Tp-z-l-N) 




(4.20) 


r(s 2 ±q- z - N) 

Here, we used the relation a-’a-^ = >S'i(j, A:)N^ with Si the Stirling numbers of first 

kind ng. The trace in ( |4.20[ ) can be evaluated using the relation 

r(l+3:±) 


r(»± - N) 


{l+x±-y±)r{y±) 

r(^±) 


k = 0 


■riSb s,{k + 1 ., + 1 ) X-V-y’ 

(4.21) 

where S 2 denotes the Stirling numbers of second kind ng and we used the abbreviations 

x± = siTp-z-l, y± = S2±q-z. (4.22) 


The relation in (4.21) can be obtained by rewriting the trace as a hypergeometric function 
which can be related to a polylogarithm using Euler’s integral transform. The series 
expansion of the polylogarithm of non-positive integer order and the integral representation 
of the Beta-function yields (4.21). The derivation can be found in Appendix From 
(4.21) it is straightforward to evaluate the sum in (4.20) using the identity 


ri+si 


*S'i(ri -f si,i)S2{i -|- 1, j -|- 1) — {j + l)5ri+sij+i + ^n+si,. 


(4.23) 


i=j 


It can be shown using the recurrence relation for the Stirling numbers of second kind and 
their inverse relation [19]. Finally, we find 


V, 


ri,r2,si,S2 


= {r2 + S 2 )! 


r(=FP=Fg + -Si--52) p±{z,q)T{Tp- Z + Si) 
r(=Fp T q + Si + r 2 + 1) P±{z,p) r(±g - z-r2) 


(4.24) 


Obviously, as Si > 0 and r 2 > 0, the normalisation in (3.7) renders a polynomial 
in z of degree si -|- r 2 . Therefore, also all entries of Q± are polynomials in the spectral 
parameter z. In the next section we will determine an npper bound for the degree of the 
polynomials in a given magnon sector. 
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4.5 Leading-order 2 ; behaviour 

In this section we show that the Q-operators are polynomials of maximal degree 
in the corresponding magnon sectors and determine the m±-dependent prefactor of the 


Q-fnnctions, cf. (2.15). 


In order to determine the leading order ^ behaviour we focus on a given intermediate 
state (ci, C 2 ,..., Cl) in the sum that yields the components of the Q-operator Q+. We 
dehne 


y{z] a, 5, c) = tr 


k+(j)(m+(j))'‘'";;“ k+(z)(m+(j)) 


b\b2...bi^ 

ClC2...Ci 


(4.26) 


cf. (4.11). Since we know the number of (z + 1 + N) in the sequenc e of a given matrix 
element of the single-row monodromies M+( 2 ;) and M+(^), cf. Section 4.3, we can expand 


ClC2...Ci 


them in powers of (;2 -|- 1). The expansion for the components of M+( 2 ;) reads 

= "V ' . (4.26) 

V J aia2...aL V J aia2...aL 

'22 ij sequences that only contain oscillators. Moreover, 


CiC 2 ...Cr 

z} ] , cf. Section 

aia2...a_t 


Here 

\ / aia2...aL 

each sequence contains the same numbers of a and a as (m+( 
and k operators N. A similar expansion holds true for the components of the monodromy 


4.3 


M+. It follows that the contribution to the sum (4.26) can be written as 


\Ar^\ 

y{z;a,b,c) = j] £ (z + 1) 

/ci=0 ^2=0 

X tr 


\4y^\+\iif'’\-ki-k2 


(4.27) 


K 


C1C2...CL 


K, 


0102 ...ai 




b\b2...b]^ 

C1C2--.CL 


When evaluating the trace using (4.24), the leading order in ^ of 3^ can be extracted 

by normal ordering the sequences in ^ ^ ^ and picking the terms with the highest 

V / 0102... 01 ^ 

My ) and picking the highest power of a. 

/ ClC2...Ci 

This is equivalent to neglecting all the quantum corrections while normal (anti-normal) 

ordering which we denote by : ^ ^ ^ : and T ^ ^ ^ T respectively. As a 

V /oi02...0i V /C 1 C 2 ...C 1 , 

consequence the large-y contribution of a given intermediate state is 


3^(y,a,&,c) = ^ ^ y '^22 

fcl=0 ^2=0 

X tr 


K^(y 


) : (mJ'^) 


ClC 2 ...Cl^ 

0102 ...Ol^ 


K 2 


OT (m 


(fc2)\^>1^2.-feL _ 
/ciC2...Ci 


(4.28) 
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where 




ClC2...C£ 

aia2...ai 


and 


M 


(fc 2 )\ 6 l 62 ... 6 L _ ^ 


+ ) 


C 1 C 2 ...CI, 



>\+k 


\4r^\+ki 


1 al"21 


^\4f^\+k2 ^\l[f^\+k2 


(4.29) 


(4.30) 


Using the trace formula in (4.24) we can evaluate the right-hand-side of (4.28). After 
neglecting further subleading contributions that arise from the trace formula we evaluate 
the sums for ki and k 2 using the relation (B.15) and hnd 

•(a,c) 


14 


(a,c)I 
ij \ 




’i+i/i;’'>i+|j>j‘>i+i4f>i 


X 


r(-p - q + 


I r(a.c)| 

I-*22 I 


+ 


1 ) 


I 7-(c,b)| 

|-'22 I 




+ 


I r(c,b) I 

|-'12 I 


21 

lA'"! 


+ \I‘ 


ic,b)\ 


21 


+ 1 ) 


T{-p - q + 


I r(“>c)| 

|-'21 I 


+ 


I r{a>c)| 

I-'22 I 


+ 


I /-(ob) I 

|-'l2 I 


+ 


I T{c,b) I 

|-'22 I 


-H) 


(4.31) 


After substituting the relation (4.16) we find that the order of z in (4.31) and consequently 
of a given magnon block of Q+ is 2m+ = + l-^i 2 ^^l + 1^2 ^|- An analogous 


derivation of the block structure holds for Q_ but it also follows from the spin-flip symmetry 
discussed in Section |4.2[ Hence, in every magnon block m± of the Q-operators Q-t the 
maximal degree of the polynomial entries is . 

Last but not least, we would like to determine the ^-independent prefactor of the 
Q-functions in (2.15). To obtain it we sum over the intermediate states (ci, C 2 ,... c^) of 
which yields the leading order in z of the Q-operator Q+ 


Q+(4 


bib2...bL ^>1 / |^(a,c)| 

aia2...aL “ 

{c}=l 


r(c,b) 


I). 


(4.32) 


Here the sum over {c} denotes the sum over all indices ci, C 2 ,..., c^. Let us introduce the no¬ 
tation riij = and dehne Cij = |/i“’^^n(/} 2 ’'^^U/ 22 ’'^^)|, where ^ {1, 2,..., L} 

simply denotes the indices whose value for the intermediate state (ci, C 2 ,..., cl) is “2”. All 
cardinalities and can be expressed in terms of riij and q,- using the relations 


HI I ~ + '^*2 ~ Cii — Ci2 , 


^(a,c)I 
H2 \ 


— Cii -)- Ci2 , 


rlc.b)! _ 


^ li 


— flu + n 2 i — C\i — C 2 i , 




(4.33) 


Hi 


— Cii + C2i 


We note that W>i in (4.31) only depends on the cardinalities and and not 

on the actual ordering of the components in the single-row monodromies. Thus we can 
parametrise the sum using the variables Cij. The sum over the intermediate states takes 
the form 


\bib2...bL ^>1 

Q + (^)j 


aia2...a£ 



3 ^(•^j Cbij , C 




(4.34) 
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From the explicit form of 3^2>i in terms of the variables introduced in (4.34) we hnd that 
the alternating sums of the binomials over Ci 2 and C 21 vanish if 77-12 7^ 0 and 7721 7 ^ 0. This 
is the case for any off-diagonal matrix element of the Q-operator. The two remaining sums 


over cii and C 22 can be evaluated using the relation (B.16). We hnd 


(q±(- 


bib2...bL ^>1 
0102...a^ 2777± 


y2m± 


’ ’ ’ S, 




^ -'01,61^02,62 ■ ■ ■ ^ 0 . 1 ^, b ^ ■ ( 4 . 35 ) 

+ P + y 

Here the leading order behaviour of the Q-operator Q_ can either be obtained from 
a similar calculation or by applying the relation among the Q-operators discussed in 
Section 14.21 

As a consequence of the results obtained in this section the eigenvalues of the Q- 
operators Q± are polynomials of degree m± in the spectral parameter z with the leading 


coefficient given in (4.35). 


5 Prom Q-operators to Q-functions 


In this section we prove that the Q-operators satisfy Baxter’s TQ-equation on the operato- 
rial level and discuss that as a consequence the roots of the eigenvalues of the Q-operators 
satisfy the Bethe equation. 


In order to show that the operators T and Q± satisfy the TQ-equation, cf. (2.14), we 


take inspiration from Baxter’s original method applied to construct the Q-operator of the 
eight-vertex model [5], see also izniEi] where a similar method has been used. The idea is 


to multiply T and Q± as given in (2.7) and (3.13) respectively and to bring the resulting 


monodromies and boundary matrices to an upper triangular form in the auxiliary space of 
the transfer matrix to hnally take the trace. To do so we use the remarkable identity 




[z = 




-1 
0+ ’ 


(5.1) 


see also j22] • Here Lq 1)_ is an upper triangular matrix in the auxiliary space of the transfer 
matrix denoted by “0” with entries in the auxiliary space of the Q-operator and the 
spin chain site i. It reads 


L, 


0+ 


.z = 


{z + lpl 

0 


1 ) 


-21 


-|- ae. 


22 


zL™(z + l) 


( 6 . 2 ) 


Furthermore, Go+ is a lower triangular matrix in the auxiliary space of the transfer matrix 
with entries in the auxiliary space of the Q-operator 


‘fO-l- 



(5.3) 


A direct computation shows that similar formulae exist for the combinations of boundary 
operators 


ICo{z)V+\2z)Ky{z) = Go+Ko+(z) 


1 — 1 
^0-h ’ 


(5.4) 
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where 


{q + z)K+{z-l) {q + z)AK+{z) 

0 ‘2z{q — z — 1 )K+( 2 : + 1) 


as well as 

where 


1-1 


^o+(^) — 

K+{z)Vl\-2z - 2)}Co{z) = Go+Ko+(z) ^o+ , 

-2{p + z){z + 1 )K+( 2 ; - 1) (j)- z -1) K+( 2 ;) a 

0 —(j) — z — l)K^{z + \ 

Using the relations above it is straightforward to show that the TQ-eqnation is satished 
die operators 

T(z)Q+(j) = 


K|,+ (2) = 


(6.6) 

(5.6) 

(6.7) 


by the operators T and Q+. From (4.6) we know that 
-1 


2z + l 
-1 

2z + l 


trotr Wo(^)L_,_( 2 ^)U+( 2 ;)K_|_( 2 ;)L+(—2^ — 2)/Co(z) 


tro tr 


Mo(z)M+(z)/Co(z)L+(2z)K+(z) 

X Mo(z)M+(z)K+(z)L+(-2z - 2)/Co(z) 


(5.8) 


where we used the explicit form of the double-row monodromies and the RTT-relation 

(6.9) 


(4.5). Furthermore, we note that 

M.(z)M+(z) = G.+ L™ 


1—1 


and 


M4z)M+(z) = G.+ LH(z)Ly-‘'(j)... V:i(z) G'l . (6.10) 


[L-1], 


[ 1 ] 


Finally, using the decompositions involving the boundary operators (5.4) and (5.6) we 
can evaluate the trace in the auxiliary space of the transfer matrix which is written as a 
product of only upper triangular matrices. We hnd 


T(z)Q+( 2 :) =tro| 




0 {q-z-l){p-z-l)^Q+{z + l) 


= {Q + z)ip + z) 


(z + l) 


2L+1 


^L+1 


-Q+( 2 : -l) + {q- z- l){p - z - 1)—^Q+( 2 ; -h 1) 


^+2 


(5.11) 


Here * denotes a non-zero entry of the upper triangular matrix. This is Baxter’s TQ- 
equation on the operatorial level involving Q+. Using the symmetry between Q+ and Q„ 
discussed in Section 4.2| and noting that 


5T(z)5-^ = T{z 


P^-P ’ 
q^-q 


( 6 . 12 ) 
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we obtain the TQ-equation for Q„. Thus we derived the TQ-equation in (2.14) for operators 


T(z)Q±(z) = {z±p){z±q) 


(z+l) 


2L+1 


2L+1 


-Q±( 2 ; -l) + {zzfp + l){z =F g + 1)—^Q±(2 + 1) . 

z + 2 

(5.13) 


From Section |4.5| we know that the eigenvalues of the Q-operators are polynomials in the 
spectral parameter of degree m± with the leading coefficient given in (4.35). Furthermore, 
we note that on the level of the eigenvalues the TQ-equation implies that 

Q±{z) Q±{-z-l) 


Q±{z-1) Q±(-z) 


(5.14) 


This is a consequence of the crossing symmetry of the transfer matrix (2.12). It follows 
that the fraction above is of the form Q±{z)/Q±{z — 1) = Yli=i^{z — zf)/{z -|- zf) with 
zt^rn± = — 1 for i = 1, 2,..., m±. Hence we obtain 


m± 


Q±(~-) = 


\[{.z-zt){,z + 4 + l). 


(6.16) 


2m± - Lzf PTi] iJi 

Finally, evaluating the TQ-equation at the roots of the polynomials Q± we hnd that the 
roots satisfy the Bethe equations (2.16). Thus we conclude that the eigenvalues of the 
Q-operators are given by Baxter Q-functions which are polynomials of degree 2m± of the 
form (5.15) with zeros at the Bethe roots. 


6 Summary and outlook 


In this article we dehned Q-operators for the open Heisenberg spin chain with diagonal 
boundary conditions. 

Our method combines Sklyanin’s formalism for integrable open spin chains |1] and 
the Q-operator construction for the closed Heisenberg spin chain HSIESIEI]. An obstacle 
that immediately arises is that the construction of the transfer matrix as presented in 
|1] assumes certain symmetries of the R-matrix and most of them are absent in the Lax 
operators used for the Q-operators construction [T3]. Thus, a priori it is not clear which 
boundary Yang-Baxter equations and unitarity relations have to be satished in order to 
build members of the commuting family of operators from the two bulk Lax matrices used 
in ra. The boundary Yang-Baxter equations and unitarity relations relevant to construct 
such operators that commute with the transfer matrix of the open Heisenberg spin chain, 
cf. Section 4T were presented in Section These relations along with the solutions to the 
boundary Yang-Baxter equations (3.6) led us to dehne two operators which commute with 


the transfer matrix. To show that the constructed operators are well dehned we studied 
the evaluation of the trace in the auxiliary space which is more involved than for the case 
of the closed chain. We found that the operators are of block structure and polynomials 
in the spectral parameter governed by the magnon number. Finally, we showed that the 
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defined operators satisfy Baxter’s TQ-equatioii on the operatorial level. This allowed us to 
deduce the form of their eigenvalues and identify them as Baxter Q-operators. 

We found that the Q-operators constructed here naturally incorporate a ^-independent 
prefactor in the corresponding Q-functions. It depends on the length, the magnon number 
and the boundary parameters p and q. In principle, this prefactor can be absorbed by 
redefining the monodromies for the Q-operator. However, naively this would spoil the 
locality which is manifest in our dehnition. 

It would be interesting to study the boundary operators introduced here in the context 
of the twisted Yangian [25] and derive further underlying functional relations for the 
Q-operators, cf. [131 Eg for the case of the closed chain. As for closed chains it should be 
possible for the open chain to obtain the Q-operators from transfer matrices by taking 
certain limits in the auxiliary space, see e.g. |2B1127] where the relevant Lax operators were 
obtained from the universal R-matrix [2H]. It may also be interesting to study how the 
Hamiltonian can be obtained from the Q-operators as studied for the closed chain in [22] 
as well as the relation to the constructions of Q-operators for the open ASEP in [3D], the 
derivation of the TQ-equation in [31] and the integral operator approach as studied for 
the closed chain in [321133] ■ 

We hope that the construction generalises to the case of non-diagonal boundary matrices. 
It was suggested in [33], see also references thereof, that Baxter’s TQ-equation is modified 
by the insertion of an inhomogeneous term. The construction of the Q-operator may shed 
more light on that relation and might yield a direct proof of the TQ-equation as presented 
here for diagonal boundary operators. 

The generalisation of the Q-operator construction to higher rank algebras is outstanding. 
In particular supersymmetric higher rank algebras with non-compact representations play 
an important role in the context of the AdS/CFT-correspondence, see [35] and references 
therein and thereof. For generalisations of the closed Heisenberg spin chain we refer the 

reader to [231 EH ED] • 

We plan to return to these questions elsewhere. 
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A Crossing symmetry of the transfer matrix 


The crossing symmetry of the transfer matrix can be shown using the crossing relation for 
the R-matrix 

0iRi2{-z-l)0^^ = -RU^), 

as well as the boundary crossing relations 

OilCi{—z — l)Oi ^ = TyyTVTT tr2 IPi2-Ri2(2^)/C2(^) , 


and 


2 (^ + 1 ) 


OiiCi{—z — l)Oi — — — tr2 Pi2-Ri 2(—2z — 2)K,2{z). 


(A.l) 

(A.2) 

(A.3) 


Here the matrix O is given by 


O = 


0 1 

-1 0 


(A.4) 


with O* = 0~^. The subscripts in the equations above denote the corresponding space. 
Using the crossing relation (A.l) and the Yang-Baxter equation (2.3) we find that (A.3) 
implies 

OiUi{-z-l)0^^ = -^^ii2^^2Ri2{-2z - 2)M\^{z)t2{z)M\^{z), (A.5) 

where lA is the double-row monodromy in ( 2.7| ) and the tq the transposition in the quantum 
space whose subscripts are suppressed. Substituting (A.2) and (A.5) into the definition of 
T{—z — 1) in (|2.7|) we find 

T(-z-l) = T(z). (A.6) 

Here we used the crossing unitarity relation in (2.5) and the relation A4q’(z) = jC{q°(z) 
with to the transposition in the auxiliary space. 


B Trace formulae 


In this appendix we evaluate the infinite and finite sums for the combination of Gamma 


functions appearing in Section 4.4 and 4.5 


In the following we repeatedly use the reflection formula of the Gamma function 


r(a;)r(l — x) = 


TT 


sin( 7 rx) ’ 

and the generalised hypergeometric function defined as 

(®l)n(® 2 )n ■ ■ ■ (®p)n Z 


0^27 • • • 7 ^p7 ^l7 ^27 ' ' ' 7 ^q7 z) ^ ^ 


where 


(c)m = c(c+1)... (c + m - 1) = 


=0 (^l)n(^2}n ■ ■ ■ (bq'jn 

r(c -I- m) 


Tic) 


denotes the Pochhammer symbol. 


(B.l) 

(B.2) 

(B.3) 
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B.l Proof of the trace formula (4.21) 


In this section we evaluate the trace cf. (4.21) which sums over the states of the infinite¬ 
dimensional oscillator space 


tr = V "i n" ■ 

r(i/± - N) ^ 


r(a:± — n) 

^0 r(i/± - n) 


(B.4) 


If /c = 0 using (B.l), (B.2) and (B.3) we can write the inhnite sum as 

^r(a:±-n) r(x±) ^ ^ 

2^ ivT—TV = Tv:rT2^i(i’ 1 - 2/±; 1 - i) • 


(B.5) 


^0 r(|/± - n) T{y±) 

The hypergeometric function 2 F 1 takes a special value if its argument z equals to 1 m 

T{c)T{c-b-a) 


2 Fi(a, 6 ;c;l) = 

r(c — a)r(c — b) 

This brings us to the result for the k = 0 case 

T{x±-N) r(l + x±) 


tr 


(B. 6 ) 


(B.7) 


r( 2 /±-N) ( 1 + a:±-|/±)r(|/±) ■ 

If k E the first term in the sum is zero, hence by shifting the summation index and 
using (B.l), (B.2) and (B.3) the sum takes the form 

E , 2,2 - !/±; 1 ,..., 1.2 - i±; 1). (B.8) 


nT()r( 2 /±-n) T{y±-1)‘ 


k-l 


Euler’s integral transform relates a generalised hypergeometric functions of a given order 
to one of a lower order via 


p+lFqJ^l (fli, . . . , Clp+1,61, . . . , feg+l, z) 


r(i,+i) 


Jo 


r(cip+i)r(&g-)-i cip-(_i) Jo 

^ ^ 2 ? • • • ; ^p-) ^ 1 ; ^ 2 ; * * * ; , 


(B.9) 


In our case we have 


fc+iFfc( 2 ,..., 2,2 |/±; 2 a;±; 1) y( 2 - L)T(v+- x^) I t) 2 ^± -± 1 


r (2 -i/±)r(j/± -x±) Jo 

X fcl 4 -i( 2 ,..., 2 ; 1 ,..., 1 ; t)dt , 


(B.IO) 


cf. (|B. 8 |). The hypergeometric function with the arguments above is related to the polylog- 

kFk-i{2 ,..., 2; 1,..., 1; t) = Li_fc(t). (B.ll) 


arithm 
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For negative integers —k the polylogarithm can be written as 


t 


U_k{t) = ^j\S2{k + l,j + 1) 
j=o 


(B.12) 


where S 2 denotes the Stirling numbers of second kind |TH]. It is the consequence of the 
repeated action of the derivative property tdtLinit) = Lin-i(f) on Lio(f) = From the 
dehnition of the Beta-function 


B{x,y)= [ — ^dt = 

Jo 


r{x)T{y) 

T{x + y) ’ 

and the reflection formula (B.l[) we obtain the hnal result 


(B.13) 


r(2/± - N) 


r(a;±) 
r( 2 /± -a;±) 


tr b.""* N* = j;(-l)P'! S2(k + 1, j + 1) 


r(2/± -x±-i-j) 


r(2/± -1 - i) 


if k e N+. 


(B.14) 


The result for the k = 0 and k G N"*" cases are exactly (4.21) 


B.2 Sums appearing in Section 4.5 


In Section [47^ two types of sums appear in the derivation of the leading order behaviour. 
The hrst reads 


^ JjY{A + k) V{A)V{B-A + j) 
t, ^ \klT{B + k) T{B-A)T{B+3)' 


(B.15) 


where A, B are parameters and n G N. For the evaluation of the sum we used (B.l), (B.2), 
(B.3), (B.6) and also the property of the Pochhammer symbol with negative integer values, 
namely that {—a)k = 0 if a G N and a < k, cf. (B.3). 

The second types of sum is 


E 


F(4-Fi) 1F(1-Fn-F4) 


^r(i + i) A r(i + n) • 

where A is arbitrary parameter and the relation can be shown by induction in n. 


(B.16) 
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